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We report a simple safe and attractive pedagogic demonstration with magnetic compasses that 
facilitates an intuitive understanding of the concept that ferromagnetism and ferroelectricity do 
not result from dipole-dipole interactions alone. Phonon dispersion relations were calculated for 
a 3-dimensional simple-cubic dipole crystal and a 2-dimensional square-lattice dipole crystal. The 
latter calculation corresponds to the compass demonstration that confirm the antiferro ground state 
structure. A sum rule for the three optical phonon frequencies is discussed. A mathematical non- 
analyticity in the longitudinal optical phonon are also illustrated. 



I. INTRODUCTION 

Ferromagnetism and Ferroelectricity in crystals have 
occupied minds of physicists and chemists for a long 
time. Both are phenomena that involve temperature de- 
pendent symmetry breaking. Ferromagnetic and ferro- 
electric states have parallel alignments of magnetic and 
electric dipoles, respectively, in all unit cells of the crys- 
tal. Other possible orderings, e.g. antiparallel alignments 
of dipoles in neighboring unit cells, correspond to anti- 
ferro states. Ferromagnetic and ferroelectric (or anti- 
ferro) states having reduced symmetries are stable only 
up to a Curie (or Neel) temperature, above which the 
crystal typically transforms to higher-symmetry para- 
magnetic or paraelectric states. Energetic stability of 
the low-temperature state reflects the resolution of com- 
petition between various interactions. Ferroelectricity 
originates from a delicate balance between short-range 
bonding and long-range dipole-dipole interactions. In 
magnets, dipole-dipole interactions are relatively weak, 
and stability is governed mostly by short-range interac- 
tions, such as direct exchange, super exchange, double 
exchange, etc. Dipole-dipole interactions in magnets are 
important, however, in the energetics of domain forma- 
tion. 

Dipole-dipole interactions vary as ^ , where r is dipole- 
dipole separation, and they are anisotropic. It is not 
straight-forward for students to intuitively understand 
the role of dipole-dipole interactions in the stability of 
macroscopically ordered states, and that ferromagnetic 
and ferroelectric states are not the results of dipole-dipole 



interactions between atoms or sites. Here, we introduce 
an attractive demonstration with magnetic compasses 
which facilitates an intuitive understanding (Sec.^. In 
Sec. mil we also introduce easy analytic and computa- 
tional calculations of the phonon dispersion relations of 
dipole crystals, which: 1) demonstrates that the dipole- 
dipole interactions yield an antiferro ground state; 2) 
yields a simple picture of LO-TO splitting; 3) reveals 
a mathematical non-analyticity at the F point in the 
phonon spectrum that arises from infinite range inter- 
actions in crystals (this last point is not readily appre- 
ciated from standard text book treatments). Finally, we 
summarize in Sec ilVI 



II. MAGNETIC COMPASS DEMONSTRATION 

If one puts two magnetic compasses on a desk well 
separated from each other, both point north as pictured 
in Fig. ^a). As one brings them closer together, their 
dipole-dipole interaction dominates and below some crit- 
ical distance, the influence of geomagnetism can be ne- 
glected. The compass needles align along the line be- 
tween them, like < SNx SN>. As in Fig. ^c), if one 
puts three, four, or more compasses in a row, one also 
finds that they align according to <SN><SN><SN> 
• • • . If, however, one arranges the compasses on a square 
lattice, they order in a stable canted antiferro structure, 
which have zero total dipole moment, as shown in Fig. [21 

We demonstrated this experiment to undergraduate 
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FIG. 1: (a) Two car windshield magnetic compasses on a 
desk, sufficiently far apart that both point north, (b) Two 
compasses sufficiently close together that the dipole-dipole 
interaction dominates, (c) Five compasses in a row. 




FIG. 2: Magnetic compasses in a 4 x 4 square lattice. 

students and had them do it themselves. It generated 
a lot of interests in them and helped them to intuitively 
understand that ferromagnetic and ferroelectric ground 
states are do not occur a result of dipole-dipole inter- 
actions alone. One author, YS, originally invented this 
demonstration^ and displayed it in the Osaka Science 
Museum. It attracted children as well as adults seeding 
interest in the origin of magnetism in crystals.^ 

This experiment is easier to perform when one uses 
magnetic compasses for car windshields, such as those 
shown in Figs. ^ and [2l Each compass has a transparent 
bubble-like plastic container filled with oil, a suction cup 
for attaching it to a windshield, and a ferrite permanent 
magnet arrow that floats freely in the oil. Dimensions of 
the ferrite magnet are 8x6x3 mm. The 8x3 mm ar- 
row sides are magnetic poles with surface magnetic fields 
^ 50 mT, that decay to 1.5 mT at the surface of the 
plastic container.^ 

I 



III. PHONON DISPERSION OF DIPOLAR 
CRYSTALS 

A system of classical dipoles on an array of sites is 
equivalent to a system of paired point charges +Z* and 
— Z* on the same sites overlappingly, in which the charges 
are allowed to move oppositely and infinitesimally: the 
dipole moment of the displaced charges is fii = 
where Ui is the displacement (separation vector) of ith 
pair point charges Along the displacement, the 

pair point charges do not shift their center of mass from 
their original site. The dipole-dipole interaction energy 
depends quadratically on the dipole moments iJ,i and is 
directly related to harmonic phonon dispersion of the sys- 
tem. Using the analysis of phonon spectrum, we show 
the stability of the antiferro structure exactly and ana- 
lytically. 

It should be mentioned that there is a difference in 
degrees of freedom between displacement Ui and a mag- 
netic compass. The number of degrees of freedom of Ui is 
three, though that of a magnetic compass is one, i.e. its 
rotation. However, as shown in the analysis of phonon 
spectrum for a 2-dimensional square-lattice dipole crys- 
tal in Sec. IIII Bl there is a correspondence between the 
structure of phonon which has the minimum energy and 
the ground state of compass lattice. 



A. 3-dimensional simple-cubic dipole crystal 

Consider a 3-dimensional simple-cubic dipole crystal 
with lattice constant ao- Dipoles are placed on lattice 
points and only the dipole-dipole interactions are consid- 
ered, 

T// A Z*'^ Ui-Uj -3{f\j ■Ui){f\j -Uj) 

V{ui,Uj;rij) = ^ , (1) 

Coo r-j 

where: Z* is the Born effective charge; eoo is the optical 
dielectric constant, refractive index squared; Ui = u{ri) 
and Uj = u{rj) are the displacements of dipoles i and 
j; Vij = Vi — Vi is di vector between dipoles i and j; and 
a hat indicates that the vector has a unit direction, e.g. 
r = r/\r\. The dipoles {u} in a super cell l^super = 
L^ao X Lyao x L^a^ = Na^ under periodic boundary 
conditions are governed by the phonon Hamiltonian: 



i=l a i=l a j=l (3 
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where ol and /3 are independent Cartesian directions (x, n is the super-ceh lattice vector: 
and z\ the a component of r runs: 



Tot = 0, ao, 2ao, • • • {L^ - l)ao 



(3) 



and ^ indicates that the summation does not include 
M* is the effective mass for dipoles, p(r,) = M*n(r,) is ^^^^^ ^j^^^j^ = ^ = q. In the 3-dimensional 
the momentum for n(r,), ^a/3(n, ) is the force constant simple-cubic lattice, direct evaluation of sum in eq. ^ is 



matrix 



difficult, because summation of in eq. (QJ converges 
slowly. However, it can be evaluated by using the Ewald 
, (4) sum technique^"^'^ as (in Z* = ecx) = 1 unit) 

I 



^a/3{rij) = [6ai3B{\rij +n|) - {vif^ n)a{rif^ n)(3Crr{\rij +n|)] 



. ^ Y^'^ikTcUk^G)^-^ exp 



G k 



\k-\-GP\ 2 2 



r 



where G is a reciprocal lattice vector, and /c is a recip- into eq. Q, the phonon Hamiltonian can be decomposed^ 
rocal vector in a first Brillouin zone of the unit cell such into the Hamiltonian which is local in /c-space. 



Ga — ' ' ' 5 —2 — , , 0, — , 2 — , • • • , (7) 

ao ao ao ao 



Lc, - 1 27r 1 27r 1 27r 1 27r 

z-Lq, ao -^q; ao -L/a (^0 ^ (^o 

(8) 

while the summation is not performed for the G = k = 
case, as indicated in eq. © by . The Ewald summa- 
tion is performed with decay functions 



Hk,a = ^^^^ + \ E ul{m.,{k)Mk), (14) 



where Ua{k)^ Pa{k), and ^a/sik) are Fourier coefficients, 
Ua{k) = ^ Ua{r) exp{-ik ■ r) , (15) 

r 

Pa{k) = ^Pair) exp{-ik • r) , (16) 



, erfcf/^r) 2 exp(— /^^r^) 

= — ^ + ^^^^s — ^ 



(9) 



and 



*a/3(fe) = Yl exp(-ifc • r) . 



_ (17) 

a,(r) = J=K(2K^ + A)exp(-«V^) , (10) 

^ ^ Here, <I>q,/3(/c) is a real function, because ^api^) is a real 

1 • • IT 1 -n. u- u J even function in the simple cubic geometry. The second 

where />: is a reciprocal decay length which one can de- jfm ^ ^ „r„ „j 

termine arbitrarily. 

It is well known that, by substituting 



= j;^^Ua{k)exp{ik-r) , (11) 



term of eq. © is transformed as 

n «0 \ 



|fc + G|^ 

4k2 



and 



(12) 



(18) 

and terms in the summation of eq. ([T^ converge rapidly 
as a function of /c + G. Therefore, with optimal n and 
fast Fourier transformation (FFT) for the first and third 
terms of eq. ^a(3{k) can be evaluated by summing 
only Tic = — L^ao, terms in real space and a G = 
term in reciprocal space. For example, if ao = 3.94 Aand 
Lx = Ly = Lz = 32, K, = 0.078 gives sufficient 
(^l^^ accuracy in ^af3{k) for double precision computations, 
as illustrated in Fig. |31 
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FIG. 3: Decay functions B{r) and Crr{r) in real space and 
exp(— A;^/K^)/ao in reciprocal space for ao — 3.94 A, Lx — 
Ly = Lz = 32, and k — 0.078 A~^. These functions decay to 
less than 10~^° at the real or reciprocal zone boundaries. 
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FIG. 4: Phonon dispersion relations for the simple-cubic 
dipole crystal along axes of high crystal symmetry in the first 
Brillouin zone calculated with a x x = 32 x 32 x 32 
super-cell. Imaginary frequencies are plotted as negative val- 
ues. 



Three eigenvalues 



< cj^ < cj? of the 3 X 3 



matrix ^a(3{k)/M^ give longitudinal-optical (LO) and 
transverse-optical (TO) phonon^^ dispersion relations of 
uj versus k as plotted in Fig. 21 where M* is the effective 
mass of a dipole. There is a sum rule, 



^1 



(19) 



because the trace of the matrix ^a/si'f^ij) in eq. Q is 
zero for all Vij and consequently the trace of ^af3{k) also 
equal to zero for all k. 

In the long wavelength limit, k 0+ (with ^ oo 
in the present super-cell calculation), ui = UJ2 converge 

This frequency 



to 



1 47rZ*2 
3 



and COS to 



47rZ*2 



M*a3 ^^^^ Y 3eooM* 

difference is called as the LO-TO splitting which has a 
non-analyticity (discontinuity) in a LO-branch at the F 
point, k = (000). The LO-TO splitting is analogous to 
the plasma-oscillation^ and corresponds to the longitudi- 
nal oscillations of dielectrics which obeys an equation of 



FIG. 5: Schematic illustration of longitudinal oscillation of a 
thin slab or film of a dielectrics. 



motion. 



M*- 



4^P. 



47rZ 



*2 



(20) 



where = —AttPz = —^TvZ'^Uz/aQ is the depolariza- 
tion field that is induced by the surface charges that are 
caused by the polarization P^, as illustrated in Fig. [51 
Along the A axis, the point group of k is C^y. Conse- 
quently, two TO phonon frequencies, coi =002^ are de- 
generate, non-degenerate ujs has to be the singly degen- 
erate LO phonon, and their ratio has to be — = ^ 
to satisfy the sum rule (|T9j) . The singly degenerate TO 
phonon at the M point has the structure of oppositely 
z-polarized rods in a checker-board arrangement on an 
xy plane. This antiferro structure has minimum energy. 
At the R point, oppositely directed dipoles align in NaCl 
structure and its three optical phonon frequencies be- 
come zero, cJi = 6^2 = ^3 = 0, because summing eq. ^ 
within n-th neighbors gives zero. For instance, within 6 
nearest neighbors, 

^ F(n,-n;ro.) = 0, (21) 

iG(6 nearest neighbors) 

and also within 12 second nearest neighbors, 

F(n,n;roi) =0 . (22) 

iG(12 second nearest neighbors) 



B. 2-dimensional square-lattice dipole crystal 

The 2-dimensional square-lattice dipole crystal, which 
has N = Lx X Ly dipoles in a fiat periodic super cell, 
is also governed by the Hamiltonian ^ and the force 
constant matrix (Q, but there are no dipoles for 7^ 
and Uz 0. In the 2-dimensional case, summation in 
eq. Q can be performed in real space. Calculated po- 
tential energies of phonons u{ri) = Uq ex.p{ik • r^) having 
wave vectors k and constant amplitudes Uq are shown in 
Fig. El Eigenmodes of polarization are Uq = (1,0,0), 
Uq = (0,1,0), and Uq = (0,0,1) for k along A axis 
and no = (1,1,0), no = (1,-1,0), and no = (0,0,1) 
for k along S axis. The TO phonon energy mini- 
mum is at the X point. There are two X points, 
k = (7r/ao,0) and k = (0,7r/ao), in the first Bril- 
louin zone of the 2-dimensional square-lattice crystal 
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FIG. 6: Energies of phonon u(ri) — u^^-x^^iih • Vj) for the 2- 
dimensional square-lattice dipole crystal, as functions of wave 
vector h are plotted, along axes M -S- F -A- X in the first 
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Brillouin zone. N 



Lx 



64 X 64 is used for this 



calculation. This dispersion corresponds to TABLE II in Ref. 
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(d) linear combination 



FIG. 7: (a) First Brillouin zone of the 2-dimensional square- 
lattice crystal. Two X points k — (7r/ao,0) and k — 
(0,7r/ao) are indicated, (b) Transverse-optical (TO) phonon 
oi k — (7r/ao,0) which is lowest in energy, i.e. u(ri) — 
(0, 1, 0) exp['i^(ri)x], is illustrated with arrows at each site. 
The direction of wavevector k is indicated with dashed arrow, 
(c) Another lowest- in- energy TO phonon oi k — (0,7r/ao), 
i.e. u(ri) — ( — 1, 0, 0) exp['i^(ri)y]. (b) A canted antiferro 
structure, which is a linear combination of (b) and (c), i.e. 
u(r,) = (0, 1, 0) exp[i^(rOx] + (-1, 0, 0) exp[i^(rOy]. 



as depicted in Fig. d^a). Corresponding two antiferro 
structures, u(ri) = (0, 1, 0) exp[z^(ri)a;] and u(ri) = 
(— 1, 0, 0) exp[i^(ri)y], are illustrated in Fig. [Tfb) and 
(c), respectively, in which oppositely polarized rows of 
dipoles are arranged alternately. A linear combination of 
these degenerate two TO phonons at different Aj-points, 
= (0,l,0)exp[i^(ri)J + (-l,0,0)exp[i^(riy, 
has a canted antiferro structure illustrated in Fig. [7fd). 
In the 2-dimentional infinite-size system, any linear com- 



FIG. 8: Comparison of energies as functions of the phase 
d of two coefficients in eq. (|23|) for finite-size 2-dimentional 
square-lattice dipole crystals from 2 x 2 to 128 x 128. 2 x 
2 configurations of dipoles are also shown in above panels 
for Q = f 5 f 5 "X' ^ (from left to right). The canted 
antiferro structure (^ = f ) gives minimum energies for all 
system sizes. Energy range becomes narrower and the energy 
converges to the value of infinite system (^ —2.549), as the 
system size gets larger. 



binations of the two lowest-in-energy TO phonons form 
a "continuously degenerate manifold of antiferro ground 
states"^. In finite-size systems, however, the continuous 
degeneracy breaks. The linear combination can be gen- 
erally expressed with two coefficients cos6> and sin^: 



u{ri) = (0,1,0) cosOexp[i — {ri)x] 

ao 

+ (-l,0,0)sin^exp[z — (rO^y] , (23) 
ao 

and, in the finite-size system, its energy becomes a func- 
tion of Lx, Ly, and 0. In Lx x Ly = 2'^ x 2^ (n = 
1, 2, • • •) system, as shown in Fig. ^0 = gives the 
simple antiferro structure (Fig. Etb)), = j gives the 
lowest-in-energy canted antiferro structure (Fig. [Zfci)), 
and = ^ the highest-in-energy structure. Therefore, 
the canted antiferro structure is realized in the square lat- 
tice of magnetic compasses as described in Sec. ^ This 
is because the canted antiferro structure minimizes leak 
flux outside the system. On the other hand, the highest- 
in-energy = ^ structure maximizes leak flux. It can 
be also seen in Fig. |S1 that, as the system size gets larger 
and consequently the finite-size effect get smaller, energy 
range becomes narrower and the energy converges to the 
value of infinite system —2.549). 

It should be additionally mentioned that, in LxX Ly = 
4x3 finite system as shown in Fig. [31 for example, there 
may be another energetics, because it has odd Ly and the 
canted antiferro structure cannot minimize leak flux out- 
side the system. Particularly in such small finite size sys- 
tems, sizes and shapes play important roles in the break 
of degeneracy. 
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FIG. 9: Magnetic compasses in a 4 x 3 square lattice. 

IV. SUMMARY 

We introduced an attractive demonstration that uses 
magnetic compasses to facilitates an intuitive under- 
standing of the fact that ferromagnetic and ferroelectric 
orderings in crystals are not the results of dipole-dipole 
interactions. It is an easy and safe experiment that may 
be a good prologue to lectures on ferromagnets and fer- 
roelectrics. Phonon dispersion of 3-dimensional simple- 
cubic and 2-dimensional square-lattice dipole crystals are 
calculated to determine which antiferro structure is most 
stable, when dipole-dipole interactions alone are consid- 
ered. For the 1-dimensional chain of dipoles, the F-mode 



with polarization along the chain is lowest in energy. For 
the 2-dimensional square-lattice dipole crystal, the an- 
tiferro X-mode with transverse polarization is lowest in 
energy. For the 3-dimensional simple cubic dipole crystal, 
the antiferro M-mode with a checker-board-arrangement 
of oppositely z-polarized rods is lowest in energy. The 
LO-TO splitting and a sum rule for the three optical 
phonon frequencies was demonstrated in a simple way. 
This is also a good example of application of the Ewald 
summation technique commonly used in simulations of 
periodic systems of charged particles. 
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